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SUMMARY 


The  epoxy  resin  Hercules  3501-6  will  be  used  as  a  matrix  material  for 
the  advanced  carbon  fiber  reinforced  structural  components  to  be  used  in  the 
F-18  Navy  fighter  aircraft. 

The  purpose  of  this  work  is  part  of  an  effort  to  determine  the  envi¬ 
ronmental  effects  (moisture  and  temperature)  on  the  properties  of  these 
composites. 

By  combining  the  experimental  and  theoretical  results  of  micro-mechanics, 
a  "finite  difference  laminate  theory,"  diffusion,  and  modeling  of  the  envi¬ 
ronment,  it  is  expected  that  a  prediction  of  the  composite  performance  as  a 
function  of  time  and  environment  can  be  made. 

The  specific  investigation  of  this  report  is  the  determination  of  the 
moisture  sorption  equilibrium  and  the  concentration  dependent  diffusion 
coefficient  of  this  resin. 

This  work  was  sponsored  by  the  Naval  Air  Systems  Command  during  FY  1978 
under  the  Task  Number  A3200000/004A/9R02200000. 


J.  R.  DIXON 
By  direction 
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INTRODUCTION 


This  work  is  part  of  an  effort  to  characterize  pertinent  properties  of 
the  Hercules  3501-6  epoxy  resin  which  is  a  material  to  be  used  as  a  matrix 
in  advanced  carbon  fiber  composites  for  structural  application  in  future 
naval  aircraft. 

Moisture  affects  the  elastic  and  strength  properties  of  organic  matrix 
composites  especially  at  elevated  temperatures.  This  is  due  to  the  fact 
that  moisture  changes  these  same  properties  in  the  resin.  It  is  therefore 
desirable  to  predict  the  state  of  moisture  sorption  with  time  and  the  ulti¬ 
mate  equilibrium  in  a  given  environment. 

The  specific  objective  of  this  work  was  to  determine  the  moisture  dif¬ 
fusion  coefficient  of  the  resin  and  the  equilibrium  concentration  as  a 
function  of  the  relative  humidity. 

Mechanical  properties  such  as  Young's  modulus,  shear  modulus,  Poisson's 
ratio,  ultimate  strain,  and  strength  have  also  been  determined  as  a  function 
of  temperature  and  moisture  concentration.  These  results  are  covered  in  a 
separate  NSWC/TR. * 


1.  Augl,  J.M.  "Moisture  Effects  on  the  Mechanical  Properties  of  Hercules 
3501-6  Epoxy  Resin,"  NSWC/WOL  TR  79-41  (in  print). 
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BACKGROUND 


A  considerable  amount  of  work  has  been  done  to  study  the  effect  of 
moisture  on  composites.  No  effort  Is  made  to  review  the  subject,  suffice 
It  here  to  refer  to  some  of  the  original  reports  and  papers. 2-13 

In  the  two  previous  reports  we  have  specifically  dealt  with  the  questions 
of  moisture  diffusion  in  composites,  by  using  a  finite  difference  approach^ 
and  with  the  question  of  how  to  model  a  natural  environment  where  temperature 
and  relative  humidity  change  constantly,  and  In  an  Irregular  fashion.  The 
results  of  these  investigations  have  led  us  to  believe  that  property  changes 
of  composites  In  a  natural  environment  under  service  or  storage  conditions 
can  be  realistically  predicted. ^ 


2.  Hertz,  J.,  et  al.  "Advanced  Composite  Application  for  Spacecraft  and 
Missiles,"  AFML-TR-7 1-186  Vol  11,  General  Dynamlcs-Convalr ,  San  Diego, 
California,  Mar  1972. 

3.  Browning,  C.E.,  and  Whitney,  J.E.  "The  Effect  of  Moisture  on  the  Prop¬ 
erties  of  High  Performance  Structural  Resins  and  Composites,"  American 
Chem.  Soc.,  Dlv.  of  Org.  Coatings  and  Plastic  Chcm. ,  33.  No.  2,  1973, 
pp.  1 37—1 48. 

4.  Augl,  J.M.,  "Environmental  Degradation  Studies  on  Carbon  Fiber  Reinforced 
Epoxies,"  Paper  presented  at  the  TTCP-Panel  P3  Meeting,  Melbourne, 
Australia,  1975,  and  paper  presented  at  the  Air  Force  Workshop  on  Dura¬ 
bility  Characteristics  of  Resin  Matrix  Composites,  Battelle-Colurabus 
Laboratories,  Ohio,  1975. 

5.  Augl,  J.M.,  "The  Effect  of  Moisture  on  Carbon  Fiber  Reinforced  Compo¬ 
sites.  II  Mechanical  Property  Changes,"  NSWC/WOL  TR  76-149,  1977. 

6.  Browning,  C.E.,  Husman,  G.E. ,  and  Whitney,  J.M.,  "Moisture  Effects  in 
Epoxy  Matrix  Composites,"  Composite  Materials:  Testing  and  Design,  ASTM 
STP  617,  1976,  p.  481. 

7.  Hofer,  K.E. ,  Rao,  N.,  and  Larson,  D. ,  "Development  of  Engineering  Data  on 
the  Mechanical  and  Physical  Properties  of  Advanced  Composite  Materials," 
AFML-TR-7 2- 20 5,  IIT  Research  Institute,  1974. 

8.  Browning,  C.E.,  and  Hartness,  J.T.,  "Effect  of  Moisture  on  the  Proper¬ 
ties  of  High  Performance  Structural  Resins  and  Composites,"  Composite 
Materials:  Testing  and  Design,  ASTM  STP  546,  1974  pp.  284-302. 

9.  Kaelble,  D.H.,  Dynes,  P.J.,  and  Crare,  L.W.,  "Interfacial  Mechanisms  of 
Moisture  Degradation  in  Graphite  Fleers,"  J.  of  Adhesion.  ]_,  25,  1977. 
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10.  Kaelble,  D.H.,  and  Dynes,  P.J.,  "Surface  Energy  Analysis  of  Treated 
Graphite  Fibers,"  J.  Adhesion.  6^,  239,  1974. 

11.  Ashbee,  K.H.G.,  and  Wyatt,  R.C.  "Water  Damage  in  Glass  Fiber/Resin 
Composites,"  Proc.  Royal  Soc..  Series  A  312.  553,  1969. 

12.  Halpin,  J.C.,  and  Pagano,  N.J.,  "Consequences  of  Environmentally  Induced 
Dilation  in  Solids,"  AFML-TR-68-395,  1969. 

13.  Pipes,  B.R. ,  Vinson,  J.R. ,  and  Chou,  T.W.,  "On  the  Hygrothermal  Response 
of  Laminated  Composite  Systems,"  J.  Comp.  Mat.,  10.  126,  1976. 

14.  Augl,  J.M. ,  and  Berger,  A.E.,  "The  Effect  of  Moisture  on  Carbon  Fiber 
Reinforced  Composites.  I  Diffusion,"  NSWC/WOL  TR  76-7,  1976. 

15.  Augl,  J.M. ,  and  Berger,  A.E.,  "The  Effect  of  Moisture  on  Carbon  Fiber 
Reinforced  Composites.  Ill  Prediction  of  Moisture  Sorption  in  a  Real 
Outdoor  Environment,"  NSWC/WOL  TR  77-13,  1977. 

16.  Augl,  J.M. ,  and  Berger,  A.E.,  "The  Effect  of  Moisture  on  Carbon  Fiber 
Reinforced  Composites.  IV  Prediction  of  Changes  in  the  Elastic  Behav¬ 
ior,"  NSWC/WOL  TR  77-61,  1977. 
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EXPERIMENTAL 


The  cured  resin  was  machined  into  plate  samples  of  .024  cm  thickness. 
The  diffusion  coefficient  and  moisture  sorption  equilibria  were  measured 
with  a  Cahn  microbalance  by  sorption  measurements.  The  details  of  the 
resin  cure  and  the  sorption  measurements  are  given  In  Appendix  A. 
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DISCUSSION 

Measurement  of  the  Diffusion  Coefficient 


The  most  direct  measurement  of  the  diffusion  coefficient  of  a  polymeric 
material  would  be  to  determine  the  steady  state  flow  of  the  diffusant  through 
a  membrane  of  the  polymer  where  the  concent  rat  ions  on  both  sides  of  the  mem¬ 
brane  are  different  but  held  constant.  The  diffusion  coefficient  is  then 
given  as  the  quotient  of  the  flow  rate  through  a  unit  area  divided  by  the 
concentration  gradient.  For  dense,  pinhole-free  films,  this  method  is  quite 
useful  and  gas  chromatographic  methods  have  been  devised  for  this  purpose. 

For  composites,  this  method  has  only  limited  value  since  it  is  quite 
difficult  to  prepare  a  thin  crack-free  composite  membrane,  though  it  may 
serve  as  a  useful  check  for  porosity,  which  is  frequently  severe  in  filament 
wound  composites.  Flow  through  pores  is  governed  by  Po isseu il le ' s  law,  i.e., 
capillary  flow  described  by  the  following  relation: 

Q  «,  a4(p1-p2)/(4ut) 

where  Q  =  volume  passing  a  circular  cross  section  of  radius  a  per  second, 

P1~P2  *  pressure  difference,  u  =  viscosity  coefficient,  i  =  length. 

Thus,  by  changing  the  pressure  difference,  it  is  easily  possible  to 
distinguish  between  diffusion  through  the  solid  polymeric  material  and 
capillary  flow. 

Other  quite  frequently  used  methods  of  determining  diffusion  coeffi¬ 
cients  of  gases  or  volatile  liquids  are  by  sorption  measurements  where  the 
sorbed  vapors  are  measured  either  volumetr ical ly  or  gravimetrically.  For 
the  experiments  described  in  this  report  we  have  used  an  electronic  micro¬ 
balance  fabricated  by  the  Ventron  Instruments  Corp.  ,  Cahn  Division. 

The  sorption  apparatus  is  schematically  shown  in  Figure  1  and  is  described 
in  Appendix  A.  The  moisture  uptake  or  desorption  is  automatically  recorded 
while  the  vapor  pressure  and  the  temperature  are  held  constant. 

In  an  earlier  report^,  we  determined  diffusion  coefficients  of  resins 
a;  ’  composites  from  the  initial  slope  of  Mf/M^  vs  t^^  (Mt  =  amount  absorbed 
at  ime  t,  M,,  *  equilibrium  concentration  at  the  respective  partial  vapor 
pressure.)  The  experiments  were  carried  out  by  exposing  rather  thick  specimens 
(5x5x0. 3  cm^)  to  some  constant  relative  humidity  and  constant  temperature 
and  weighing  the  samples  from  time  to  time.  This  method  is  simple,  however 
time  consuming,  and  therefore  not  very  useful  if  absorption  and  desorption 
measurements  are  to  be  carried  out  until  final  equilibrium  is  reached. 
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FIGURE  1  MOISTURE  SORPTION  BALANCE 
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In  our  present  investigation  we  were  interested  in  whether  or  not  the 
3501-6  epoxy  resin  had  a  concentration  dependent  diffusion  coefficient.  It 
was  therefore  desirable  to  obtain  the  total  rather  than  just  the  initial 
sorption  curves. 

The  concentration  dependence  of  the  diffusion  coefficient  is  usually 
obtained  by  performing  sorption  or  permeation  measurements  at  different 
boundary  conditions  on  geometrically  well-defined  samples  for  which  the 
diffusion  equation 


3£ 

3t 


(where  c  =  concentration,  t  =  time,  and  D(c)  is  the  diffusion  coefficient 
which  may  depend  on  the  local  concentration  of  the  diffusant)  can  be  solved. 
Closed  form  solutions  for  Fick's  equation  are  available  only  for  concentra¬ 
tion  indpendent  diffusion  coefficients.  If  D  depends  on  concentration, 
Fick's  equation  has  to  be  solved  by  numerical  methods. 

Since  it  may  become  important  to  compare  experimental  data  among  vari¬ 
ous  laboratories,  it  might  be  helpful  here  to  discuss  briefly  the  essential 
aspects  of  determining  the  diffusion  coefficient  from  sorption  measurements 
on  plate  samples.  For  a  sheet  of  infinite  extension  and  exposure  to  the 
same  vapor  pressure  on  both  sides,  the  following  closed  solution  has  been 
derived  for  a  diffusion  with  concentration  independent  D. ^ 


M„/M  =  1  -  l  - r  exp[-D(2n+l)2A/4Jl2] 

„-0  <2n+1>  ’ 


•2(?y 


'  71/2+2 


ierf  c 


and  for  "longtime"  exposure 


VM« 

t  co 


- exp  [-D(2n+l)2TT2t/  f  2] 

(2n+ir 


where  Mt  and  M«,  are  the  amount  of  moisture  taken  up  after  a  time  t  and 

at  equilibrium  (i.e.,  after  infinite  time),  t,  =  sheet  thickness.  Formula 

(2)  can  be  used  to  derive  the  useful  formulas  (4  and  5)  for  determining 

or  D?  or  D. since  the  summation  term  in  (2)  vanishes 
l  l  1/2 


17.  Crank,  T. ,  The  Mathematics  of  Diffusion.  2nd  ed. ,  Oxford  Univ.  Press, 
1975.  ~~~~ 
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where  R 


?  1/2 

d(Mt/Mj/d(t/r) 

D1/2  -  . 0/»9/(t/*2)  1/2 


(A) 


(5) 


Che  subscript  i  indicates  that  l)  was  obtained  from  the  initial  portion  of  the 
sorption  curve.  The  subscript  1/2  indicates  the  time  after  1/2  of  the  equi¬ 
librium  concentration  has  been  absorbed  or  desorbed.  In  other  words,  if  the 
diffusion  coefficient  does  not  depend  on  concentration  the  results  obtained 
from  absorption  or  desorption  experiments  are  the  same.  If  the  diffusion 
coefficient  depends,  however,  on  concentration  then  D,  obtained  from  such  a 
measurement,  is  really  some  average  d i f fus ion_coe f f ic ient  which  we  have,  in 
view  of  t^e  following,  indicated  with  a  bar  (D).  Also,  then  the  equality 
between  D  and  D  no  longer  holds. 

From  equation  (3)  another  formula  (6)  can  be  derived  which  is  useful 
for  "longtime"  experiments,  when  the  summation  can  be  effectively  replaced 
by  its  first  term  from  which  (6)  and  (6a)  follow: 


M 

In (i- 


Dll 2 1 


(6) 


[d  «.n(Ma)-Mt)]dt  -  ir25/i2  ^ 

From  this  it  follows  that  a  logarithmic  plot  of  (1-M^/bO  vs  time  becomes 
1  inear . 

In  some  cases  it  may  be  impractical  to  wait  for  equilibrium  to  be 
attained.  For  such  cases  another  useful  formula  (7)  may  be  derived  in  a 
similar  way  from  (3). 

a  (  r  5n2M  1 

<Mt>2  "  (Vi  "  ^2  exp!.~  e2  J 


-  (Mt>1  exp 


(•¥*]! 


(7) 


12 


NSWC  TR  79-39 


The  concentration  dependent  D(c)  can  be  obtained  by  solving 

C 

<Co>  -  <!>  V5'3/  °  c2/3,HC,dC  (8) 

J  o 

as  described  by  Crank^  or  by  solving 

rCo 

Dn(C  )  =  1.85(C  )-1*85  /  (C  -C)0,85  D(C)dC  (9) 

o 

1 8 

as  described  by  Kishimoto  and  Enda  .  Both  methods  require  a  series  of 
absorption  and  desorption  curves  respectively.  Graphical  or  numerical  dif¬ 
ferentiation  is  carried  out  to  give  D(c). 

19 

Frensdorff  approached  the  solution  of  Kick's  equation  (where  D(c)  is 
concentration  dependent)  by  a  numerical  method. 

He  has  shown  that  the  slopes  of  the  linear  portion  of  curves  of  equa¬ 
tion  (6a)  where  D  is  replaced  by  D(c )  are  identical,  i.e.,  -n^/4.  Some 
authors  have  called  the  average  diffusion  coefficient  obtained  from  the  short 
time  solution  as  "integral  diffusion  coefficient"  while  that  obtained  from 
the  long  term  solution  as  "differential  diffusion  coefficient."  These  names 
indicate  that  the  short  term  solutions  were  obtained  in  a  process  where  the 
diffusion  coefficient  has  a  wide  range  due  to  the  steep  concentration  gradient 
while  the  long  term  solution  requires  a  measurement  close  to  the  equilibrium 
where  the  gradient  is  very  flat  and  therefore  the  range  of  D(c)  is  small. 

We  have  therefore  used  the  Frensdorff  method  for  our  measurement  wh ich  makes 
a  numerical  integration  according  to  equations  (8)  or  (9)  unnecessary  and 
D(c)  is  obtained  directly.  D(o),  the  diffusion  coefficient  as  zero  concen¬ 
tration  is  approached,  is  also  obtained  directly  from  a  desorption  experiment 
rather  than  by  an  iterative  method  as  described  by  Crank^O. 

Duda  and  Vrentas^l  have  described  a  complicated  method  (method  of 
moments)  to  obtain  the  total  range  of  D(c)  from  a  single  sorption  curve.  The 
authors  have  indicated,  however,  that  highly  accurate  longtime  sorption  curves 
are  essential  for  this  method.  The  authors  did  not  compare  the  results  of 
their  theoretical  model  with  actual  experimental  data. 


18.  Kishimoto,  A.,  and  Enda,  Y. ,  J.  Polym.  Set  A.  1799,  1963. 

19.  Frensdorff,  H.K. ,  J.  Polym.  Set.  A.  2 ,  341,  1964. 

20.  Crank,  T. ,  and  Park,  G.S.,  Trans.  Faraday  Soc.,  43.  240,  1949. 

21.  Duda,  T.L. ,  and  Vrentas,  T.S.,  A1CHE  Journal.  Mar  1971,  p.  464. 
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Errora  Due  to  Experimental  Inaccuracies 


1.  Effect  of  Temperature  Control 


Over  a 
with  temperature 


amall  temperature  range  the  diffusion  coefficient 
as  indicated  by  equation  (10) 


changes 


D  ■  Dq  exp 


(10) 


a 8  a  rule  of  thumb,  D  increases  by  about  a  factor  of  2  for  every  10*C.  This 
means  that  an  error  in  the  temperature  measurement  of  1*C  gives  an  approxi¬ 
mate  6Z  error  in  the  measured  diffusion  coefficient.  If  an  average  tempera¬ 
ture  of  a  small  temperature  fluctuation  (+1*C)  can  be  ascertained,  the  error 
should  be  considerably  smaller  (less  than  1Z).  There  is  a  small  error  intro¬ 
duced  if  the  experiment  is  carried  out  on  plate  samples  that  have  to  be 
removed  froa  the  exposure  chamber  to  be  weighted.  The  error  is  largest  at 
the  initial  t  is**  periods  but  gets  successively  smaller  at  the  longtime  expo¬ 
sure  because  of  the  square  root  dependence  of  moisture  uptake  with  time. 
Therefore,  a  continuous  measurement  of  moisture  uptake  in  a  closed  system 
by  means  of  a  quartz  spring  balance  or  a  Cahn  electro  balance  is  more  desir¬ 
able  if  accurate  data  are  required.  On  the  other  hand,  the  sample  size  for 
such  instruments  is  quite  small  (100-300  mg). 


2.  Effect  of  Sample  Thickness 


The  thicknesses  of  the  sheet  samples  range  usually  from  .250  to  .025 
cm  for  composites.  If  a  micrometer  of  +.001  cm  accuracy  is  used  for  the  thick¬ 
ness  measurement,  then  the  errors  in  thickness  measurements  range  from  .A  to  A 
percent  and  the  corresponding  diffusion  coefficients  range  in  error  from  .16  to 
16  percent  since  the  thickness  enters  as  a  square  in  equations  (A)  to  (7). 

3.  Effect  of  Relative  Vapor  Pressure  and  Weight  Measurements 


The  error  caused  by  inaccuracies  in  vapor  pressure  measurement 
can  be  twofold: 


a.  The  error  in  D(c)  is  proportional  to  its  slope, 

b.  The  error  is  also  dependent  on  the  accuracy  with  which  the  pres¬ 
sure  can  be  held  constant  during  the  final  phase  of  the  sorption  or  desorption 
experiment  for  determining  the  slope  that  enters  equation  (6a).  Similarly 
important  is  the  accuracy  of  the  weight  measurement  which  must  be  better  than 
1Z  for  the  long  term  solution.  The  combined  error  in  D  may  reach  M 2Z . 

A.  Effect  of  Time  Measurement 

Since  most  sorption  experiments  require  times  from  several  hours  to 
many  days,  time  should  be  the  least  source  of  error  (assuming  the  recorder  has 
been  calibrated). 
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5.  Range  of  Usefulness  of  Equation*  (4)  to  (7)  for  Determining  the 

1)  if  fun  ion  Coefficient 

The  equations  (4)  to  (7)  apply  strictly  only  for  determining  diffusion 
coefficients  that  do  not  depend  on  concentration.  Even  then  the  range  is 
limited  because  only  one  or  two  terms  of  the  infinite  series  equation  1-3  were 
used.  Formula  (4)  gives  less  than  It  error  in  the  range  where  Mt/M.  is  less 
than  0.5.  The  same  accuracy  is  obtained  with  formula  (5).  Equation  (b)  gives 
less  than  1  percent  error  in  the  range  Mt /M„  >  0.7.  In  using  formula  (7), 
the  error  depends  on  both  Mj  and  and  on  the  difference.  For  instance,  if 
Mi  is  such  that  Mt/M«  "0.5  and  M2  such  that  Mt/M„  -  .97,  the  error  is  less 
than  1.6  percent.  If  Mj  and  M2  are  both  chosen  to  be  such  that  Mt/M.  >  0.8 
the  error  is  less  than  1  percent. 

6.  Edge  Effects  of  Plates  with  Finite  Dimensions 

A  thorough  analysis  of  the  edge  effects  for  various  geometries  was 
given  by  Rothwell  and  Marshall^, 

The  correct  ion  formula  for  the  short  time  solution  of  isotropic 
materials  is  * 


D-D 

o 


i+i+£h 

«  b  J 


(where  D  -  corrected  diffusion  coefficient,  Dc  -  diffusivitv  for  a  plate 
with  infinite  dimensions  (obtained  from  equation  (4)),  1  -  plate  thickness, 
a  "  length,  b  -  width).  For  the  long  term  solution  (equation  (6)  and  (7)) 
Kothwell  and  Marshall  suggested 


2.  -,-1 


“  -  no  [U2  (V  j) 

L  it  K  J 


for  a  and  b  >>  1 ,  and  Mt /M„  >  0.6.  These  equat ions  were  applied  for  the 
i)A,  5^,  D*(c)  and  D^(o)  given  in  this  report. 


22.  Rothwell,  W.S.,  and  Marshall,  H.P.,  "Analysts  of  Experimental  Transport 
Data  Diffusion  of  Water  In  EPDM,"  Lockheed  Missile  4  Space  Co.  Technical 
Report,  LMSC -0566642,  Oct  1977. 

23.  Shen,  C.H.,  and  Springer.  G.S.,  .1.  Comp.  Mat..  10.  2,  1976. 
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RESULTS 


The  specimen  used  for  the  sorption  and  desorption  experiments  had  the 
following  dimensions:  2.5  x  1.8  x  0.024  cm-*. 

The  sample  was  exposed  in  the  sorption  apparatus  at  the  desired  tem¬ 
perature  and  pressure  and  the  moisture  gain  or  loss  was  recorded.  A  typical 
absorption  (A)  and  desorption  (D)  curve  given  by  a  short  term  plot  (equa¬ 
tion  (4))  is  shown  in  Figure  2.  The  same  data  plotted  according  to  equa¬ 
tion  (6a)  are  shown  in  Figure  3.  From  the  straight  line  portions  of  these 
plots  were  determined  the  6A,  69,  D^(c)  and  D^(o),  values,  listed  in  Table  1. 

From  the  difference  in  the  sorption  and  desorption  curves  it  is  quite 
apparent  that  there  must  be  a  concentration  dependent  diffusion  coefficient. 

Figure  4  shows  a  plot  of  Mt/M.  vs.  (D0t/h2)*/2.  if  the  diffusion 
coefficients  were  constant  and  equal  to  2.08E-10  cm2/sec  as  found  for  D*?(o) 
at  50*C,  the  absorption  and  desorption  curve  should  be  identical  and  follow 
the  curve  indicated:  (theoretical).  Using  the  same  scale,  one  observes  that 
the  experimental  curves  differ  considerably. 

Figure  5  shows  a  graph  of  the  resin  diffusion  coefficient  at  30*  and 
50*C  as  a  function  of  moisture  concentration  in  weight  percent.  (From  a 
resin  density  of  1.274  this  can  be  readily  converted  in  volume  percent  if 
so  desired. ) 

Figure  6  shows  the  ratio  of  D(c)/D(o). 

The  equilibrium  concentration  of  moisture  in  the  resin  with  respect  to 
the  surrounding  relative  humidity  is  shown  in  Figure  7.  For  an  AS-carbon 
fiber  composite  of  a  fiber  volume  fraction  of  0.62  the  calculated  equilibrium 
concentration  of  moisture  is  indicated  by  the  broken  line.  An  actual  sample 
exposed  at  80Z  RH  absorbed  1.1539  percent  of  water  which  agrees  quite  well 
with  the  calculated  moisture  uptake  of  1.187  (for  Vf  «  0.62,  composite 
density  ■  1.5825,  and  resin  density  ■  1.274). 

Figure  8  shows  the  same  data  in  terms  of  moisture  activity  in  the  resin 
versus  volume  fraction.  Since  at  chemical  equilibrium  the  activity  aj  of 
moisture  is  equal  to  the  ratio  of  p/p0  (where  p  is  the  partial  vapor  pres¬ 
sure  and  p0  is  the  saturation  pressure  of  moisture  over  liquid  water  at  the 
given  temperature)  there  is  a  simple  relationship  of  moisture  activity  at 
equilibrium  and  relative  humidity  at  which  the  experiment  was  performed: 

ZRH  ■  100  x  aj.  (Such  a  plot  is  also  called  a  sorption  isotherm.) 
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FIGURE  3  MOISTURE  ABSORPTION  ANO  DESORPTION  OF  HERCULES  3501-6  AT  50°  C  AND  80%  RH  (L  -  .024  cm). 


Table  1.  Diffusion  Coefficients  of  Hercules  3501-6  Epoxy  Resin  (Corrected  for  Edge  Effects) 
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sec  was  taken  as  an 


SORPTION  CURVE 


FIGURE  4  THEORETICAL  AND  EXPERIMENTAL  SORPTION  AND  DESORPTION  CURVES 
FOR  3601-6  RESIN  (DIO)  -  2.08  E  -10)  cm2/**  AT  50°  C  AND  80%  R.H. 
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MOISTURE  CONCENTRATION.  (WT  %) 

FIGURE  6  CONCENTRATION  OF  DEPENDENCE  OF  THE  MOISTURE  DIFFUSION  COEFFICIENT 
OF  HERCULES  3501  6  EPOXY  RESIN. 


NSWC  TR  79-39 


23 


NSWC  TR  79-39 


From  Figure  5,  it  appears  that  the  moisture  diffusion  coefficient  goes 
through  a  maximum  at  least  at  ambient  temperature.  Such  a  behavior  is  not 
unknown  and  it  may  occur  with  diffusants  that  can  readily  form  hydrogen 
bridges  and  therefore  lead  to  self  condensation. 

If  such  a  self  condensation  (clustering)  of  a  solvent  in  a  polymer 
occurs,  the  diffusion  of  the  solvent  will  be  reduced  since  the  size  of  the 
diffusant  has  increased  to  form  a  dimer  or  even  larger  condensation  clusters. 
The  answer  of  whether  or  not  such  a  condensation  occurs  may  be  obtained  from 
statistical  thermodynamics  of  polymers. 

Based  on  the  investigation  of  McMillan  and  Mayer2^  on  the  statistical 
thermodynamics  of  multicomponent  systems,  Zimin  has  proposed  two  forms  of 
"clustering  funct ions . "25-27 


This  is  to  be  understood  that  the  mean  number  of  type  1  molecules  (=H20)  in 
the  neighborhood  of  a  given  type  1  molecule  in  excess  of  the  mean  concentra¬ 
tion  of  type  1  molecule  is  $iGn/Vi  ■  CjGji,  where  and  Cj  are 

the  volume  fraction,  partial  molar  volume  and  molar  concentration  of  compo¬ 
nent  1,  respectively. 

This  clustering  function  is  easily  derived  from  the  sorption  isotherm 
(Figure  8),  from  which  aj$i  vs  a^  is  obtained  (Figure  9  and  Table  2).  A 
graphical  differentiation  of  the  curve  of  Figure  9  gives  finally  G^j/Vj  or  C^G]^ 
The  latter  function  is  shown  in  Figure  10.  As  a  comparison,  the  clustering 
function  of  unstretched  Nylon  66  is  also  shown,  where  the  sorption  data  were 
taken  from  the  literature.2® 


24.  McMillan,  W.G.,  and  Mayer,  J.,  J.  Chem.  Phys..  13.  276,  1945. 

25.  Zimm,  B.H. ,  J.  Chem.  Phys..  14.  164,  1945. 

26.  Zimm,  B.H.,  ibid.  21,  934,  1953. 

27.  Zimm,  B.H. ,  and  Lundberg,  J.L.,  J,  Phys.  Chem..  52.  425,  1956. 

28.  Bull,  H.B.,  J.  Amer.  Chem.  Soc..  66,  1499,  1949. 
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RELATIVE  VAPOR  PRESSURE,  (a,) 


FIGURE  9  «/atVERSUS  a, OF  HERCULES  3601-6  EPOXY  RESIN 


Table  2.  Zimin's  Clustering  Functions  Gu/Vj  and  CiGu  of  Hercules  3501-6  Epoxy  Resin 
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CLUSTER  FUNCTION. 
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Since  CiC}}  is  the  mean  number  of  water  molecules  near  a  given  water 
molecule  in  excess  of  the  mean  concentration,  the  mean  size  of  the  cluster  is 
CiGii+1,  which  indicates  that  the  size  of  the  cluster  increases  rapidly 
at  higher  relative  humidity  equilibrium  concentrations. 

Similar  results  have  been  reported  in  the  literature  by  Starkweather^ 
for  cellulose,  keratin,  collagen,  serum  albumin,  polyvinyl  acetate  and  poly¬ 
methyl  methacrylate. 


29.  Starkweather,  H.W.,  Polymer  Letters.  _1,  133,  1963. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


1.  The  moisture  diffusion  coefficient  of  Hercules  3501-6  has  been  deter¬ 
mined  at  30'  and  50'C  as  a  function  of  concentration. 

2.  It  was  found  that,  at  very  low  moisture  concentration,  the  diffusion 
coefficient  is  low,  while  at  concentrations  between  1  and  5X  of  moisture  (equi¬ 
valent  to  30-90Z  relative  humidity  equilibrium  concentration),  the  diffusion 
coefficient  is  rather  flat  and  seems  to  have  a  maximum  at  lower  temperatures 
(30'C).  This  behavior  can  be  rationalized  from  a  clustering  mechanism  of  water 
at  higher  moisture  concentrations.  Such  clustering  reduces  the  dlffusant 
mobility  and  thus  its  diffusion  coefficient. 

3.  The  concentration  dependence  was  directly  determined  from  the  long 
term  sorption  experiments  according  to  Frensdorff's  method. 

4.  It  is  recommended  that  diffusion  studies  also  be  made  on  stressed 
resin  and  composite  samples  in  order  to  obtain  information  on  whether  or  not 
these  materials  show  stress  Induced  diffusion. 
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APPENDIX  A 


Sample  Preparation 

The  resin  Hercules  3S01-6  was  obtained  from  the  manufacturer  (Hercules 
Bacchus,  Utah)  and  kept  stored  In  a  freezer. 

The  resin  was  heated  to  93'C  at  which  temperature  it  was  degassed  under 
vacuum.  Then  the  liquid  resin  was  cast  into  a  plate  mold  and  cured  under  the 
following  conditions:  the  resin  was  kept  for  1  hour  each  at  the  following 
temperatures:  100',  115',  124',  130',  and  150'C.  The  final  cure  was  then  10 

hours  at  177'C.  The  resin  was  then  allowed  to  cool  and  was  machined  into  thin 
(0.024  cm)  plate  samples. 


Sorption  Chamber  and  Sorption  Measurements 

Figure  1  shows  a  schematic  of  the  sorption  apparatus  which  consists  of 
a  Cahn  electrobalance,  a  weighing  chamber,  a  vapor  generator,  a  capacitance 
manometer  with  a  range  from  0.1  to  1,000  mm  of  Hg  pressure,  a  constant  tem¬ 
perature  bath,  a  5  buffer  vapor  reservoir,  a  vacuum  pump,  and  a  recorder  with 
a  sensitivity  of  0.002  mg.  The  temperature  in  the  wooden  Isolated  chamber  was 
usually  held  at  3'C  higher  than  in  the  sample  chamber  by  air  circulation  with 
two  blowers.  The  sample  chamber  was  held  at  constant  temperature  by  circulat¬ 
ing  water  through  a  heating  jacket. 

The  predried  sample  was  transferred  into  the  weighing  chamber  and  vacuum 
was  applied  with  the  valves  V},  V^,  and  V3  open.  The  water  temperature  of 
the  constant  temperature  bath  was  raised  to  80'C  to  remove  moisture  that  may 
have  been  sorbed  during  the  sample  transfer.  Then  the  temperature  was  lowered 
to  the  desired  measurement  temperature.  Valve  V2  was  closed  and  valve  V4  was 
opened  till  the  desired  vapor  pressure  was  reached  (as  indicated  by  the  pres¬ 
sure  gauge).  Instead  of  adjusting  the  vapor  pressure  directly  from  a  reservoir 
of  liquid  water,  it  is  also  possible  to  keep  a  binary  salt  water  mixture  in  the 
vapor  generator  which  restricts  the  relative  pressure  to  only  certain  fixed 
values  (dependent  on  the  kind  of  salt  used  and  on  the  measurement  temperature). 
The  sorption  was  run  till  equilibrium  was  attained. 

At  relative  humidities  higher  than  80Z  the  Cahn  balance  became  unstable. 
Whenever  possible  it  is  desirable  to  do  the  entire  series  of  sorption  and 
desorption  measurements  without  removal  of  the  sample.  It  took  about  100 
seconds  till  the  pressure  was  reached  at  which  the  measurement  was  to  be 
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carried  out,  especially  at  the  high  relative  vapor  pressures.  (If  by  acci¬ 
dent  the  desired  measurement  pressure  Is  exceeded,  this  can  be  easily  recti¬ 
fied  by  opening  valve  V2  to  the  vacuum  pump.)  It  Is  also  desirable  to  be 
able  to  control  the  valves  from  outside  the  insulation  box. 
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